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Abstract. We review aspects of black hole thermodynamics, and show how entanglement 
of a quantum field between the inside and outside of a horizon can account for the area- 
proportionality of black hole entropy, provided the field is in its ground state. We show that 
the result continues to hold for Coherent States and Squeezed States, while for Excited States, 
the entropy scales as a power of area less than unity. We also identify location of the degrees 
of freedom which give rise to the above entropy. 



1. Introduction 

The field of black hole thermodynamics started around the time Bekenstein and Hawking 
argued that black holes have entropy and temperature, given respectively by [1]: 



where M and Ah are the mass and horizon area of the black hole, and ipi = \J Ghjc? is 
the Planck length in four dimensions. Together, they satisfy the following laws, which are 
analogous to the laws of classical thermodynamics ^ : 

d{M(?) = TudS + + VLdJ (First law) (2) 
^Sbh > (Second law) . (3) 

The question naturally arises as to whether there are miscroscopic, elementary states, consistent 
with the black hole temperature and other macroscopic parameters, whose logarithm gives its 
entropy: 

Sbh = \^^. (4) 



^ In subsequent formulae, we will assume c = 1 = /i. 



and if so, then what those states are. Another related question in this area of research is 
the following: is the area proportionality of the entropy, also know as the Area Law (AL) in 
Eq.(l) (as opposed to volume proportionality, which is the case for normal thermodynamic 
systems) more generic, or in fact universal, beyond the realm of black holes? This has given 
rise to the so-called Holographic Hypothesis, which conjectures that the information content of 
our universe may be encoded in a two dimensional screen which surrounds three dimensional 
space. 

Finally, the expectation that black holes, after absorbing matter (and information contained 
therein) will radiate itself away in the form of pure thermal radiation at the temperature (1), has 
given rise to the Information Loss Problem: thermal radiation does not contain any information, 
except that of the temperature of the emitting source, so if the black hole disappears completely, 
does the information do so as well? While there is nothing in principle which prevents this form 
happening, it requires the quantum mechanics which governs the evolution to be non-unitary, 
since unitary evolution takes pure states to pure states. This eventually raises the question: 
does quantum mechanics break down in the presence of a black hole? To be able to address 
one or more of the above questions, one has to first explore origins of Black Hole Entropy 
and temperature. Among the various approaches that have been undertaken, the two most 
popular are: String Theory and Loop Quantum Gravity. In the first, degeneracy of strings 
attached to D-branes are claimed to give rise to the required degeneracy, while in the second, 
the spin-network states at the horizon are invoked for the same purpose. There have been some 
other interesting approaches as well [2, 3]. 

In addition to the above, there has been another idea, which goes as follows: consider a free 
scalar field, propagating in a spacctimc containing a black hole. Now, do a partial tracing over 
the degrees of freedom (DOF) of the field which are inside the horizon. The resultant density 
matrix is mixed, due to the information being 'hidden' by the horizon (although the entire 
system, before tracing, is pure). When one computes the Entanglement Entropy associated 
with this, the latter turns out to be proportional to the horizon area [4, 5]! [The analytical 
proof for the area law have recently been shown in Refs. [6].] It is this approach which we 
will follow in this article. Here we will critically examine an important assumption that the 
authors had made in their analyses: that the scalar field is in its Ground State (GS). While this 
simplifying assumption is justified and even natural to make in a first calculation, we would like 
to see how important it really is to validity of the AL. Indeed we will find that while replacing 
the GS by generic Coherent State (CS) or a class of Squeezed States (SS) does not affect the 
AL, doing so with a class of Excited States (ES) does change the AL, and in fact renders the 
entropy proportional to a power of area less than unity. 

One can ask as to why we consider a scalar field. While we do not have a completely 
satisfactory answer to this question, we are motivated by the fact that a significant part of 
gravitational pertTirbations in black hole backgrounds are scalar in nature [7]. So the scalar 
field does not have to be introduced by hand from outside. 

This article is organised as follows: in the next section, we briefly review entanglement in 
quantum mechanics and its role in entropy. In Section 3, we set up the formalism for computing 
entropy of entangled scalar fields in flat spacetime, when the field is in its ground state. In 
Section 4, we extend our results to generic coherent states and a class of squeezed states. In 
Section 5, we examine excited states and their effects on entanglement entropy. In Section 6, 
we try to pin-point the location of the DOF which are responsible for entanglement entropy. 
In Section 7, we extend our results to black hole spacetimes. In Section 8, we conclude with a 
summary and open questions. 



2. Entanglement entropy 

Consider a quantum mechanical system, which can be decomposed into two subsystems u and 
V, as shown below ^: 



System = Subsystem u + Subsystem v 



Correspondingly the Hilbert space of the system is a Kronecker product of the subsystem 
Hilbert spaces: 

H = nu®ny . (5) 

If \ui) and \vj) are eigenbases in Hu and Hv respectively, then \ui) ® \vj) form an eigenbasis in 
H, in terms of which a generic wavefunction l^') in H can be expanded: 

l"^) = ^dij\ui)0\vj) en . (6) 

ij 

Note however, that in general, |^) cannot be factorised into two wavefunctions, one in u and 
the other in v. 

I*) 7^ (8 . (7) 

Such states are called Entangled States or EPR states. The ones that can be factorised, on 
the other hand, are called Unentangled States. For example, for a system with two spins, the 
following is an unentangled state: 

I Ti) + I TT) = I T) ® (li) + I T)) , (8) 

while the following is an entangled state: 

IU) + ITT)7^I---)®I---) ■ (9) 

Entangled states have a variety of uses, including in Quantum Teleportation. 

Next, let us define density matrices. If the quantum-mechanical wavefunction of a system 
is known (however complex the system might be), its density matrix is defined as: 

P=l*)(*l- (10) 
It can easily be verified that p satisfies the following properties: 

IpI > , p^ = p, p^ = p. (11) 

The last property is known as idempotency, from which it follows that the eigenvalue Pn of the 
density matrix can only be or 1. Thus the Entanglement Entropy or Von Neumann Entropy, 
defined as follows, vanishes : 

S = -Tr{p\np) = -'^Pnlnpn = . (12) 

n 

Now, one can take the trace of p, only in the subsystem v, to find the Reduced Density Matrix, 
which is still an operator in subsystem u: 

Pu = Try{p) = ^{vi\p\vi) = dijdlj\ui){uk\ . (13) 

I i,k,j 
^ for details, we refer the reader to the review [8]. 



For p„, it can be shown that the following properties hold: 

\pu\ > , pi = pu , pI^Pu ■ (14) 

That is, idcmpotcncy no longer holds. As a result, its eigenvalues now satisfy: < < 1, 
and the entanglement entropy is non-zero: 

Su = -Tvu {Pu In Pu) = -^ Pn(u) lnp.„(„) > . (15) 

n 

The ignorance resulting from tracing over one part of the system manifests itself as entropy. 
One important property of reduced density matrices is that if we traced over u instead and 
found the reduced density matrix p^, then the latter would have the same set of non-zero 
eigenvalues. Consequently, the entanglement entropies are equal, being a common property of 
the entangled system: Sy = Su- 



3. Entanglement entropy of scalar fields 

Now consider a free scalar field in (3-|-l)-dimensional flat space, with the following Hamiltonian: 



7r^(x) + \V(p{x 



Decomposing the field and its conjugate momentum into partial waves 



Im 



Im 



yields: 



1 r°° { 

Im Im ^ K 



dr 



(16) 



(17) 



(18) 



Next, discretise along the radial direction with lattice spacing a, such that r — > = 
ia; rj+i — ri = a. The lattice is terminated at a large but finite N ('size of the universe'), 
and an intermediate point n is chosen, which we call the horizon, and which separates the 
lattice points between an 'inside' and an 'outside'. As of now, there is nothing intrinsically 
special about this point, although later we will associate it with a black hole horizon. The 
discretisation is depicted below: 



n(Horizon) 



N 



Inside Outside 
The discretised version of Hamiltonian (18) takes the form: 



Im 



2a 



N 



7.2 I (i I ("PlrnJ VlmJ+l V 1(1 + I) 2 



(19) 



Note that it resembles the Hamiltonian of a set of N coupled harmonic oscillators, with 
the interaction between them contained in the off-diagonal elements of the matrix Kij (the 
coordinates Xi replace the field variables ifim)'- 



1 ^ 1 ^ 



(20) 



The density Matrix, tracing over the first n ol N oscillators {x = Xn+i, ■ ■ ■ , xn), is given by: 

/n 
JJ dXii;{xi,...,Xn;Xn+l,...,XN) ll}*ixi,...,Xn;Xn+l,...,x'pf) (21) 
i=l 

where {x. = Ux , UKU^ = Diagonal). The wavefunction 



N 



^p{xi, ...,xn) =YlNi H^^ik^) -xA exp I ■ Xj 



1=1 



(22) 



denotes the most general eigenstate of Hamiltonian (20), being a product of N harmonic 
oscillator wavefunctions. Note that ^ pout, since pout is mixed, although the full state 
is pure. Consequently, the entanglement entropy = — Tropin p) > . For the scalar field 
Hamiltonian (19), K^j can be read-off: 



1] 



l{l + l)6ij + - dadji +[N--] 6iNSjN + P + o + r 



(23) 



Note that the last two terms denote (nearest-neighbour) interaction and originates in the 
derivative term in (18). Schematically, 
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(24) 



where the off-diagonal terms represent interactions. The GS of the above Hamiltonian is given 
by: 



N 



i){xi,...,XN) = JjA^iexp(-^A;|).xf) 



(25) 



which follows from (22), by setting all Vi to zero. The corresponding density matrix (21) can 
be evaluated exactly (the suffix signifies GS): 



Pq{x, x') ~ exp [— (x"^7x + x''^ ^x')j2 -|- x^ (ix\ 



(26) 



where: 



A B 
C 



(27) 



Note that B and /? are non-zero if and only if there are interactions. The gaussian nature of 

the above density matrix lends itself to a series of diagonalisations ^ , such that it reduces to a 
product of {N — n), 2-oscillator density matrices, in each of which one oscillator is traced over 
[5]: 



N-n 



+ PiViv'i 



The corresponding entropy is given by: 

N-n 



s 



5] -ln[l-e,]-^lne. 
i=l ^ ^* 

Thus, for the full Hamiltonian H = J2im ^im, the entropy is: 

oo 

s = Y,m + i)Si 



Pi 



1 + 



(28) 



(29) 



(30) 



where the degeneracy factor (21 + 1) follows from spherical symmetry of the Hamiltonian. In 
practice, we will replace the upper bound of the sum in the above to a large value Imax- For 
the interaction matrix (23), the above entropy, computed numerically, turned out to be [4, 5]: 

2 



5 = 0.3(n + l/2)^ = 0.3{ - 



(31) 



The logarithm of the entropy versus the logarithm of the radius of the horizon is plotted in 
Figure 1, and the slope is seen to be 2. 



4. Coherent and Squeezed States 

4-1. Coherent state 

In this section, we first examine the situation when GS in (25) is replaced by coherent states. 
A CS ('shifted' ground state) is defined as the following, for which ApAx = h/2: 

JL 1 1 

iIj{xi,...,xn)cs = [[Niexp{--kl)i{^i- aif) (32) 

i=l 

= [[exp{-i^^ai)Niexp{~-kfy.x^) (33) 
i=i ^ 

where j)j = —idjd^^ . Next, defining x = x — U~^a , dx = dx, it is easy to show that the 
density Matrix for the CS can be written as: 

/n 
Yl dxi 'il^csixi,...;xn+i,...) i^hsi^i^-'-'^^n+i,---) = p{x;x') . (34) 

i=l 

3 V-rV^ = diag, = iD^VpV'^-t'K WPW^ = diag, Vi € v = {V-iV'^)^VT . 



log(S) = 1 .9904 log(n + 0.5) - 1 .2097 
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log(R) 



Figure 1. In(S') vs ln(i?) for GS. The slope is 2. 



Remarkably, the density matrix has the same function form as the GS density matrix (26), albeit 
in terms of the tilded variables. However, this means that the eigenfunctions will be identical 
to those for the GS, as a function of the tilded variables. Consequently, the eigenvalues, and 
hence the entropy will be the same as those of the GS [10]! 

4-2. Squeezed state 

SS are also minimum uncertainty packets, for which for which ApAx = h/2. However, they 
are characterised by Ap ^ 1 or Ax <C 1 (or vice- versa). Here, we consider a class of SS, 
characterised by a single parameter r, and given by the following wavefunction: 



N 



ipssixi, . . .,xn) = r^/^ JJiVj exp 



i=l 



1/2 2 



(35) 



In this case, defining a new variable as: x = -y/r x , dx = ^/r dx. , the density matrix becomes: 



Pss [x;x 



/n 
dxiii)ss{xi,---\Xn+i,...) ^ssi^i^-'-'^^n+i,---) = p{x\x') . (36) 
i=l 



Once again, it has the same functional form, and the entropy is same as that for the ground 
state [10]. 



5. First excited state 

Although, ideally one would like to compute the entanglement entropy for the most general 
eigenstate (or superpositions of such eigenstates), given by Eq.(22), this turned out to be a 
formidable task. Instead, to examine the role of excitations which do not form minimum 
uncertainty packets, we consider the following wavefunction, which is a superposition of N 
wavefunctions, each of which signifies exactly one oscillator in the first excited state ^. 



N 



ip^{xi . ..Xn) 



J2 aiNiHi k 



1 



9 D] -. 



(37) 



* The analysis for two harmonic osciUators was done in [9]. 



= V2 [a^Klxj ^o{xu...,XN,) (38) 

Where = (ai, . . . , ajv), and normalisation requires a^a = 1. The density matrix can be 
written as: 

/n 
JJ dxi [x'^Ax] -00 {xf, x) ■00 {xf, x') . (39) 
i=i 

This can be evaluated exactly, yielding: 

1 



1 - - {x^Ajt + x'^Ajx') + x^Apx' 



Poix,x') (40) 



tT tA „ „r 7^3 TT - f Ab 



p{x,x') 

where 

A = [/^ Klaa' KlU^^l^^ J (41) 

Without further approximations, (40) cannot be factorised into 2-particle density matrices. 
However, noting that po in the above is a gaussian (given by (26)), decaying virtually to zero 
beyond its few sigma limits, we check whether quantities ei,2, defined below, are small for any 
given A matrix. That is, whether: 

^1 — XjYKxx^pXmax ^ 1 > ^2 = X^y^^^A^X^ax ^ 1 (^^) 

where 

corresponding to 3(7 limits in the argument of po- If both €1,2 <C 1, then we can make the 
following approximation: 



1 — - {x^ A-yX + x'^^A-yx') + x'^ Ajjx' exp 



-- [x'^A^x + x'^A^x') + x^ A^x' 



(44) 



with which (40) can be written as: 

p(a;, x') K exp ^- ^ (x^Yx + x'^Yx') + x^/^'x'^ [/?' = /? + A^, 7' = 7 + A^] . (45) 

The above being of the same form as the gaussian density matrix, but with shifted parameters, 
can be factorised once again into two particle density matrices, and the associated entanglement 
entropy can be evaluated. The rest of the computation is numerical (done using MATLAB). 
We choose the following values for parameters: N = 300, n = 100 — 200, o = 10 — 50, where 
o signifies the number of non-zero entries in the vector: = (1/y^) (0, • • • , 0, 1, ■ • ■ , 1 ). For 
the above choice, it was seen that ei 2 < 10-3 . We give the results of our computation in the 
form of relevant graphs. Figure 2 shows the logarithm of the entropy for the GS and ES for 
o = 10, 30, 40, 50. We see that for the ES, the entropy scales as with lesser a for higher o. 
It is less than unity for any o > 0. The AL docs not seem to hold! Next, in Figure 3, wc plot 
(2/ + 1)5"; vs / for the GS as well as for various ES and several n. It is seen that for the GS, 
there is a peak at ^ = (s-wave), followed by another one around / 50 — 100, the degeneracy 
factor of {21 + 1) being responsible for the latter. For ES, the first peak shifts to Z > 0, with the 
shift being greater for larger o. The graphs do not show a second peak, although there seems 
to be an increase towards higher values of I. Thus, higher partial waves are seen to get excited 
with greater excitations. 



ln(S )Vsln(n + 1/2) for = 10 




log(S) = 1 .8957 log(n + 0.5) - 0.7458 
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log(S) = 1 .651 log(n + 0.5) + 0.61 58 
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log(S) = 1.7695 log(n + 0.5) - 0.0337 
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log{S) = 1.5218 log(n + 0.5) + 1.3498 
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Figure 2. Plot of ln(5) vs ln{R) for ES. The slope is less than 2. 
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Figure 3. Plot of {21 + 1)5/ vs I for GS and different ES for various n {N = 300). 



6. Where are the degrees of freedom? 

In this section, we address the question: where are the DOE, which give rise to entanglement 
entropy for black holes?^ 

To this end, we re-consider the interaction matrix K in (24), and set to zero, by hand, the 
interaction (off-diagonal) elements except for a small window (of say size 3x3), and move the 
window from the origin through the horizon, to the outer edge of the spherical box, as shown 
below: 

^ In Ref. [11], the authors have studied the scaUng properties of energy fluctuations close to the boundary. 



K 



I ^ 
I X 



X 
X 
X 



(46) 



x| 
x| 



V 



With the above modified matrix, now we compute the entanglement entropy, and plot it as a 
function of the position of the centre of the window (the green dotted boxes). The result is 
shown in Figure 4 (the solid red curve). Note that when the window lies entirely outside or 
inside the hrizon, the entropy is zero. When the window just enters the horizon, the entropy 
starts to rise, and when it is symmetrically placed between the inside and outside, the entropy 
rises to 100% of its value, and decreases again as it is further moved. The above results confirm 
the intuition that entanglement between inside and outside the horizon is the source of this 
entropy, and shows how it rises from zero to its maximum value as the entangled DOF are 
gradually incorporated. In Figure 5, we probe this further, by keeping one extremity of the 
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Figure 4. Entanglement entropy comes from the near-horizon DOF. 



window just outside the horizon, while increasing its width from zero to the horizon radius. 
We see that with a width of one lattice spacing, about 85% of the total entropy is obtained, 
and within a width of about 8 — 10 lattice spacings, it reaches 100%. While this shows that 
most of the entropy comes from DOF very close to the horizon, a small part (about 15% in 
this case), has its origin deeper inside. The latter is understood to indirectly affect the DOF 
near the boundary via the nearest-neighbour terms in (24). We expect a similar analysis for 
ES to give us a better understanding of the scaling of the entropy with area in that case. 
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Figure 5. Rise of entropy as interaction region is increased. 



7. Generalisation to black hole spacetimes 

In this section, we generahse the framework of our calculations, so that they are applicable to 
black hole spacetimes. We start with the Schwarzschild metric with a horizon at r = tq: 



f{r) r 



and define the following coordinate transformation from {t, r) to (r, R) coordinates: 
T = t + ro 



In 



1 - \A7n) 



1 + \A7n) 



2J — 



2 ri 



2 

3 1 

^3 



' 



such that (47) transforms to the following metric, in Lemaitre coordinates: 



ds' 



-dT^ + 



dR^ 



+ 



2ro 



4 

3 2 



The Hamiltonian for a free scalar field in the above spacetime can be written as ^: 



1 f°° 
H{t) = - I dR 



(47) 



(48) 



(49) 



3{R-t) 2 

Next, choose a fixed Lemaitre time, say r = and perform the following field redefinitions: 

Mr) 



(50) 



vr(r) = ^/riTi{r) , (j){r) 



(51) 



the current analysis is done for I = 0. Generalisation to any Z > is straightforward. 



Then, at that fixed time, the Hamiltonian (50) transforms to: 



1 r 



dr 



(52) 



Note that this is identical to the Hamiltonian in flat spacctime, Eq.(18). Also, it follows from 
Eq.(48) that the horizon now corresponds to R = |ro. Thus, one is free to trace over either 
the region R = ^ |ro or over the region R = |ro —>■ oo, and all our results go through for 
a fixed r. Once again, the AL holds for GS, CS and SS, while it does not for ES. We also 
show that although the degrees of freedom near the horizon contribute most to the entropy, 
the DOF deep inside must also be taken into account for the AL. 



8. Summary and open questions 

In this article, we have shown that entanglement between DOF inside and outside the horizon 
is a viable source of black hole entropy. However, while the AL holds for the GS, CS, and SS, 
this entropy scales as an area power less than unity, when the considered wavefunction is a 
superposition of first ES. Furthermore, while the DOF near the horizon contribute most to the 
entropy, those farther away give a small contribution as well. We would like to investigate this 
further for ES, as that might provide a physical understanding of the deviation from the AL. 
We would also like understand the phhenomena of Hawking radiation and information loss in 
this icture. We hope to report on this elsewhere [12]. 
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